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The dynamical behaviors of FRW Universe containing a posivive/negative potential scalar field 
in loop quantum cosmology scenario are discussed. The method of the phase-plane analysis is used 
to investigate the stability of the Universe. It is found that the stability properties in this situation 
are quite different from the classical cosmology case. For a positive potential scalar field coupled 
with a barotropic fluid, the cosmological autonomous system has five fixed points and one of them 
is stable if the adiabatic index 7 satisfies < 7 < 2. This leads to the fact that the universe just 
have one bounce point instead of the singularity which lies in the quantum dominated area and 
it is caused by the quantum geometry effect. There are four fixed points if one considers a scalar 
field with a negative potential, but none of them is stable. Therefore, the universe has two kinds of 
bounce points, one is caused by the quantum geometry effect and the other is caused by the negative 
potential, the Universe may enter a classical re-collapse after the quantum bounce. This hints that 
the spatially flat FRW Universe containing a negative potential scalar field is cyclic. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The cyclic model of universe was first noticed by Ein- 
stein who considered the possibility of this model as an al- 
ternative to the model of an expanding universe. But this 
model is failed for the entropy problem jl]]. Since then, 
many cyclic universe models had been introduced, i.e., 
the Brane cosmology model 0, [3, the Baum-Frampton 
model [4|, cyclic universe caused by loop quantum grav- 
ity (LQG) [a-Q' 6tc. The universe may enter a cyclic 
stage if it contains a scalar field with a negative poten- 
tial [lol . [T]| . Considering the flat Friedmann-Robertson- 
Walker(FRW) universe, the Friedmann equation reads 



quantum gravity (LQG) [14l.[l5|. In the present isotropic 
and homogeneous settings, the canonically conjugate 
fields of L QG are parameterized by constants c and p 
respectively p^: 



c I, 



-1 0, 



and El ^ p Iq ' 



(2) 



(1) 



with K = 8ttG, and (j) is the scalar field. Obviously, if 
the potential V{(j)) is negative, the universe may enter 
a re-collapse phase, i.e., the Hubble parameter H — 0, 
and then a = if = -¥{((>) with 1/(0) < 0. But 
there are other possibilities that the universe will not 
enter the cyclic phase because the quantum singularity 
can not be avoided sometimes. To avoid the quantum 
singularity, we consider the universe described by loop 
quantum cosmology (LQG). 

LQC [13, [13 is a canonical quantization of homoge- 
neous spacetime based on the techniques used in loop 
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It has been shown that the loop quantum effects can 
be very well described by an effective modified Fried- 
mann dynamics. Two corrections of the effective LQG 
are always considered: the inverse volume correction and 
the holonomy correction. But the inverse volume modi- 
fication suffers from gauge dependence which cannot be 
cured and thus yields unphysical effects. Moreover, for 
a Universe with a large scale factor, the inverse volume 
modification to the Friedmann equation can be neglected 
and only the holonomy modification is important. 

In this paper, we will use phase-plane analysis of 
the cosmological autonomous system to investigate the 
dynamical behavior of the FRW universe with a posi- 
tive/negative potential scalar field (see the Eq. ([4])) in 

effective LQG based only on the holonomy modification. 
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By this, the Friedmann equation adds a — term to 
the right-hand side of the standard Friedmann equation 
p^ . Since this correction comes with a negative sign, 
the Hubble parameter H, and then a will vanish when 
p = Pc, and the quantum bounce occurs. Based on the 
holonomy modification, the dynamical behavior of dark 
energy has recently been investigated by many authors 
(17.-.22i] . The attractor behavior of scalar field in LQG 
has also been studied [23l - l25l |. The dynamical behavior 
of scalar field with positive/negative potential(Eq. ^) 
in classical cosmology has been discussed by [lOl • 

The organization of this paper is as follows. In the Sec. 
im we present the general equations of motion (EOM) . In 
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the Sec. IIIIl we introduce the autonomous system and 
analyse the dynamical behavior of this system. The nu- 
merical results will be presented in the Sec. IIVI Finally, 
the conclusions will be devoted in the Sec. |Vl 



II. BASIC EQUATIONS 

We focus on the flat FRW cosmology. The modified 
Friedmann equation in the effective LQC with holonomy 
correction can be written as llal 



1 



p_ 

Pc 



(3) 



in which p is the total energy density and the natural unit 
K = 8ttG = 1 is adopted for simplicity. We consider a 
self- interacting scalar field (j> with a potential Vicj)) (for a 
more general consideration of self-interacting scalar field 
theory in the universe with holonomy effects resulting 
from the Loop Quantum Cosmology, see [l^) coupled 
with a barotropic fiuid. Then the total energy density 
can be written as p = + pj, with the energy density 
of scalar field p^ = ^(fP' -I- V{4)) and the energy density 
of barotropic fiuid p^. It is worth pointing out that for 
the case of negative potential if the energy density of 
the scalar field remains positive, must satisfying 0^ > 

-V{4>)- 

We consider a simple exponential potential for the 
scalar field 6 



-X4> 



(4) 



where A is a dimensionless constant characterizing the 
slope of the potential ^ld\ , and Vq can be a positive or a 
negative constant corresponding to positive or negative 
potential. The positive/negateive potential has been dis- 
cussed by many authors in classical cosmology [iq | and in 
loop quantum cosmology. If Vq > 0, the energy density 
of scalar field p is always positive, and the Hubble con- 
stant H will be zero if and only if p = pc, which is just 
the quantum bounce caused by the quantum geometry 
effect. But if Vb < 0, it is obvious that the Hubble pa- 
rameter H can be zero if p = 0, i.e., pj + ^(fP = —¥{(/)). 
This means that the universe has one quantum bounce 
phases when p — p^, and one re-collapsing phase lying in 
the classical stage when pj + = —V{(l)). 

We consider that the energy momentum of this field to 
be covariantly conserved. Then one has 



'(j) + 3H<j) + V' = 0, 
p^ + i^Hp^ — 0, 



(5) 
(6) 



where 7 is an adiabatic index that satisfies p-y = (7 — l)p^ 
with being the pressure of the barotropic fiuid, and 
the prime denotes the differentiation w.r.t. the field (j). 
Differentiating Eq. ([3]), and using Eqs. ^ and ([6]), one 
can obtain 



H = - 



+IP1 



1 - 



2(P7 + P-p) 



Pc 



(7) 



Considering the situation of Vo < 0, obviously, H > 
when p = Pc, while H <0 when p^ + = —V{(j)). Lets 
assume the universe is expanding at the beginning with- 
out losing generality, then the scalar factor a is increas- 
ing and the total energy density p is decreasing. When 
p = 0, i.e., + P7 = ~V{(I)), = and < 0, the 
turnaround happens. The universe ceases expanding and 
turns to a contracting phase. In the contracting phase, 
the scale factor a is decreasing, while p is increasing. In 
this phase, if p — pc, H = 0, H > 0, the bounce oc- 
curs, and the universe reenters the expanding phase. To 
conclude, a scalar field with negative potential in loop 
quantum cosmology can evolve cyclically. We will show 
these pictures in the Sec lIVI Noted that the analysis in 
this paper is based on effective LQC. This is different 
with [9,], which is based on the quantum theory of LQC. 



III. DYNAMICAL ANALYSIS 

In the last section, we introduced the basic equations. 
From Eq. we can see that a new term is added to 
the Friedmann equation in effective LQC. Based on the 
research in [2^, this lead to the consequence that the 
dynamical system will be very different from the one in 
classical cosmology i-^-' O'^^ additional dimension 

will appear in the autonomous system of LQC. So, the 
stability properties of the universe will change. It is in- 
teresting to see if there exists any stable point in the new 
autonomous system and checking out whether the result 
is consistent with our prediction presented in the Sec. 
im because if the autonomous system contain a stable 
point, the universe will finally evolve into a steady state 
and will not cyclic. This discussion also can help us to to 
understand how LQC changes the conclusions obtained 
in classical situation. In this and the next section, we will 
discuss the dynamical behavior of the universe containing 
a scalar field with positive/negative potential. 

Equations ©-([S]) characterize a closed system which 
can determine the cosmic behavior. We introduce three 
dimensionless variables ^1 



^ _ VWl _ p 

— ;= — 1 V — ^= — , z— — 
V6H VSH Pc 



(8) 



Then the Friedmann equation (see Eq.Q) can take the 
simple form 



and one can get 



_Pj_ 



1- z 



2 2 



(9) 



(10) 



We use upper/lower signs to denote the two distinct cases 
of ±V > here. The z term is a special variable in LQC. 
< z < 1 when the Universe lies in quantum dominated 
stage, i.e., p < pc, and z = when the universe lies in 
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classical stage, i.e., p <ti Pc- Differentiating Eq. ([9]) one 
can obtain 



barotropic fluid in the next two subsections and show 
some numerical analysis in the next section. 



H 



3.^ + ^4 



1 



1- z 



2 2 

X TV 



(l-2z). (11) 



Based on the above dimensionless variables ([8]), we 
will use the phase-plane analysis of the cosmological au- 
tonomous system to investigate the dynamical behavior 
of this exponential potential scalar field coupled with a 



A. Autonomous system 

Using the new variables ([8|), and considering Eq. (fTO 
and Eq. ([TT|) . we can rewrite Eqs. ([3|)-([5|) in the flowin; 
form 



— — — —3a; -f- X 



dx 
dN 

dy 
dN 

dz 
dN 



X (l-2z)±^/|A^/^ 



.2 . 37. 1 



3x' + ^{ 



= 2z(l-z) 



2 '1 - z 



-3x^ 



X (1 - 2z) - \I-Xxy, 



2 '1 - z 



x'Ty') 



(12) 
(13) 
(14) 



where iV = In < 



Obviously, the terms on the right-hand side of Eqs. 
(fT2)) - (fT4)) only depend on x,y,z but not on N or other 
variables. Such a dynamical system is usually called an 
autonomous system. We will restrict our discussion of 
the existence and stability of critical points of this au- 
tonomous system to the upper half-plane y > 0, consider 
only A > without loss of generality. 

This autonomous system is 3-dimensional for we 
choose the exponential potential (|4]). If one do not con- 
sider this special potential, it is possible to get one or 
more extra dimensions [23j [ . 



According to this condition, it is easy to get the values 
of the fixed points, and they are shown in the Table I. 

The properties of each fixed point are determined by 
the eigenvalues of the Jacobi matrix 



M 



/di df df\ 

dx dy dz 

dg dg dg 

dx dy c 

dh dh dh 

dx dy dz , 



(16) 



B. Stability properties 



Now we consider the fixed points of this autonomous 
system described by Eqs. ([11])- (HH). Let / = dx/dN,g = 
dy/dN, h — dz/dN . The fixed points of the autonomous 
system satisfy the condition 



0. 



(15) 



According to Lyapunovs linearization method, the sta- 
bility of a linearized system is determined by the eigenval- 
ues of the matrix M . If all of the eigenvalues are strictly 
in the left-half complex plane, then the autonomous sys- 
tem is stable. If at least one eigenvalue is strictly in the 
right-half complex plane, then the system is unstable. 



Remarks on Tab. I: 

(i). In fact, for y < 0, we have another fixed point 
however, for simplicity, we just consider < 7 < 2. This 



means that yc is a complex number and is thus a non- 
physical solution. Since the adiabatic index 7 satisfies 
< 7 < 2, all the terms that contain 7 should not change 
sign. If one consider 7 = or 7 = 2, some of the eigen- 
values corresponding to points Pi, 2, 3, 5 will be zero. To 
analyze the stability of such points with zero eigenvalues. 
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TABLE I: The stability analysis of an autonomous system in LQC. The system is described by a self-interacting scalar field 
(j) with exponential potential V coupled with a barotropic fluid p^. Explanation of the symbols used in this table: Pi denotes 
the fixed points located in the 3-dimensional phase space, which are marked by the coordinates {xc, yc, Zc, Ac). M'^ means the 
inverted matrix of the eigenvalues of the fixed points. U stands for unstable, and S stands for stable. 



Fixed-points 




yc 


Zc 


Existence 


Eigenvalues 


Stability 


Pi 











AU V,X and 7 


M'^' = (-3 + §7, -37, §7) 


U, for all A, 7 


P2 


1 








AU V,X and 7 


M'^ = {-6, 6 - 37, 3 - ^A) 


U, for all A, 7 


Ps 


-1 


1 





AU V,X and 7 


M'^ = (-6,6 - 37,3-h ^A) 


U, for all A, 7 




A 

V6 







V > and A^ < 6 


MT = (-A^-3+4,A^-37) 


S, for A^ < 37 


P4_ 


A 

V6 







1/ < and A^ < 6 


M^^ (a^-^,-3+^,A^-37) 


U, for all A, 7 




2A 


\/^(2-7) 





> and A^ > 37 


See Eq.fTTl) 


U, for all A, 7 



we need to resort to other more complex methods [23 
(ii). The eigenvalues of the Point P5 is 



tive/negative potential by using numerical tool. 



-37 



M = 



3 -2A+7A+-y/4A2-20A27+9A272_2473+4872 



3 -2X+i\-^4X^ -2OA27+9A272 -247^+487^ 
^4 A 



(17) 



Obviously, 4A^ - 
be hel d. This 
7A+y4A2 

-2A + 7A 



- 20 A27 9 A272 - 2473 + 487^ > should 
point is unstable because -Fi = — 2A -I- 
4872 < and F2 = 

_ 2473 + 4872 < 



2OA27 



9A272 



2473 



^4A2 - 2OA27 + 9A272 
can not be satisfied at the same time. It is also impossible 
that Fi = F2 = 0, since this means that 7 = 2. We do 
not consider this condition. 

Based on Table I and the remarks above, we can find 
that there is no stable point when one consider the scalar 
field with a negative exponential potential. There is just 
one stable point (the Point P4+) when A2 < 87 and 
A2 < 6, if one consider the scalar field with a positive 
exponential potential. This conclusion is different from 
the stability properties of the scalar field with positive 
exponential potential in classical cosmology [27j . 

We find the values of coordinates (x^ yc) of P4^ and 
the fourth fixed point in the Table I of [27| , which is in the 
classical setting, are the same. So are their stability prop- 
erties. We also can find there exists one more stale point, 
i.e., the fifth fixed point in the Table I of It has the 
same values of fixed point with P5, but obviously, in our 
case, P5 is unstable. This difference comes from the mod- 
ification of Friedmann equation. Comparing the classical 
Friedmann equation Eq.([T|) and the modified Friedmann 
equation Eq. (jS]) , we can find that the modified Friedmann 

equation adds a — 3p term, and then the dimension of 
the autonomous system in LQC is different from the one 
in classical cosmology. 



IV. NUMERICAL RESULTS 

In this section, we will analyze the dynamical behav- 
ior of FRW universe containing a scalar field with posi- 



y 0.6 




FIG. 1: The phase picture of the positive potential projection 
on the X — y plane. The barotropic fiuid is radiation, i.e., 
7=1- 

We found that there is not any stable fixed points if 
one considers a scalar field with a negative exponential 
potential, and there is only one stable fixed point if one 
discuss a scalar field with a positive exponential poten- 
tial. We presented this stable fixed point in FigH] As it 
shows, almost all solutions of the dynamical equations of 
the autonomous system end up at P4_^ . 

The evolution pictures of H{t),a{t), 4>{t), p{t), (j){t) and 
P'y{t) are shown in the Figl^J We consider the time 
t = is the quantum bounce point. Then one has 
H(t = 0) = 0. The initial conditions of /C7 and a 
are (p{t = 0) = l,p^i{t = 0) = 0.01 and a{t = 0) = 0.1 for 
both the universe containing a scalar field with positive 
potential and the one with negative potential. According 
to these initial conditions and the Friedmann equation 
(see the Eq. (|3])), one can get the initial condition of (j). 
It is easy to find that the spatially fiat FRW universe 
is cyclic if we consider the scalar field with a negative 
exponential potential, but not cyclic if one consider the 
positive exponential potential. We have to point out that 
the initial conditions do not affect the cyclic behavior of 
the universe, different initial condition only change the 
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locations of the fixed points and the universe evolution 
rate. 

Here, we analyze these evolution pictures focusing on 
the situation of negative exponential potential. Fig. 2(a)- 
(c) show that the evolution of the Hubble parameter 
H{t), the scale factor a{t) and the total energy density 
p. Just like we pointed out before, there are two suf- 
ficient conditions to insure H{t) can be zero. The one 
takes place in the contracting phase, where the scale fac- 
tor is decreasing, total energy density is increasing, and 
when p — pc, H — 0, the bounce occurs. Then the 
universe enters into an expanding phase. Different from 
the classical cyclic universe, with loop quantum correc- 
tions, the universe bounces back before a singularity oc- 
curs. The other one occurs when p = 0, which means 
i</)^ -I- V{4>) + Pj = 0. Because the total energy den- 
sity equals to zero, this condition lies in the classical 
stage. From the image, we see this corresponds to the 



zero points where H(t) evolves from positive to nega- 
tive, a(t) evolves from increasing to decreasing, and p{t) 
evolves from decreasing to increasing. From Fig. 2(d) and 
Fig. 2(e), we can see that the scalar field is not cyclic, and 
its momentum continuously decreases. This consequence 
comes from the modified Friedmann equation. Combined 
with Fig. 2(c) and Fig. 2(f), we can find that (j){t) can not 
reach exactly zero, so that (j){t) will not turn around but 
will continue to increase and just slow down during the 
bounce and re-collapse. The bounce and re-collapse in 
scale factor and Hubble parameter are not the bounce 
and re-collapse in 4>{t). So there is a series of bounces 
and re-collapses at finite values of </>(t). Meanwhile, since 
the potential is not constant, the behavior of bounce and 
re-collapse of the universe is not periodic but still cyclic 
^5]. Just as Fig. 2(a) shows, the cycle period of H{t) is 
rising gradually. 



H(t) 0.0 




FIG. 2: The evolution pictures of the universe containing a positive potential/negative potential scalar field, for representing 
(a) H{t), (b) a{t), (c) p{t), (d) cj){t), (e) 4>{t) and (f) p^{t), respectively. The initial conditions are H{t = 0) = 0, a{t = 0) = 0.1, 
pj{t = 0) = 0.01, (j>{t = 0) = i, and (p{t = 0) = 1.4922 for positive potential (the pointed line), and H{t = 0) = 0, 
p-f{t = 0) = 0.01, (j}{t = 0) = |, and ^(i = 0) = 1.007 for negative potential (the solid line). The barotropic fluid is radiation, 
i.e., 7= |. 



V. CONCLUSIONS ponential potential scalar field coupled with a barotropic 

The goal of this paper is to discuss the dynamical prop- 
erties of the FRW universe with a positive/negative ex- 
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fluid in LQC scenario. 

The Friedmann equation is modified in LQC scenario. 
In this situation the universe will enter a quantum bounce 
phase when the total energy density p is equal to the 
critical energy density pc- After quantum bounce, the 
universe will immediately enter a super-inflation phase, 
and then enter a stage which is described very well by 
general relativity. According to the phase-plane analy- 
sis of the evolution of the spatially flat FRW universe 
in LQC scenario, we found that there is just one stable 
fixed point (see the Point in the Table I). So, if the 
universe is described by flat FRW metric, just includ- 
ing the scalar field with positive potential and coupled 
with a barotropic fluid, it has only one quantum bounce. 
This is different from the stability properties of classical 
cosmology which includes the same fields |27|| . The differ- 
ence comes from the modification of Friedmann equation. 
Since the Friedmann equation is modified in LQC sce- 
nario, the autonomous system is 3-dimensional in LQC, 
while 2-dimensional in classical cosmology. 

Also, we have presented a phase-plane analysis of the 
evolution of a spatially flat FRW universe containing a 
scalar field with a negative exponential potential V = 
Vb exp(— A(/)) (Vb < 0) plus a barotropic fluid. Because of 
the different potential, the dynamical equation and thus 
the dynamical behavior of these two situations are differ- 
ent. When we consider the situation of the scalar field 
with positive exponential potential, the autonomous sys- 
tem has five fixed points, and one of them is stable if 



one constrains the adiabatic index with < 7 < 2. But 
when the potential is negative, there exist just four fixed 
points, and none of them is stable under the same 7. The 
universe may enter a classical re-collapse after the quan- 
tum bounce in this scenario if -f p-y = — Vb exp(— A(/)). 
The spatially flat FRW Universe is cyclic. 

To summarize, from the mathematical point of view, 
the necessary condition for the universe become cyclic is 
that there is no stable point. In the model we discussed, 
just because the positive potential case does have one sta- 
ble point, almost all solutions of the dynamical equations 
of the autonomous system end up at the stable point. 
This means finally, the universe will evolve into a very 
steady state, rather than a cyclic stage. Meanwhile the 
negative exponential potential is not the only potential 
that can lead to cyclic universe. For a cyclic universe 
model caused by another kind of scalar field potential in 
LQC, see [t^, while 81 considered the situation where the 
space is anisotropic. 
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